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Abstract 



We study an integro-differential equation that describes the slow erosion of granular 
flow. The equation is a first order non-linear conservation law where the flux function 
includes an integral term. We show that there exist unique traveling wave solutions that 
connect profiles with equilibrium slope at ±00. Such traveling waves take very different 
forms from those in standard conservation laws. Furthermore, we prove that the traveling 
wave profiles are locally stable, i.e., solutions with monotone initial data approaches the 

' traveling waves asymptotically as t — > +00. 

> \ 

keywords: traveling waves, existence and stability, integro-differential equation, conser- 
vation law. 

> ■ 1 Introduction 

On 

■ We consider the Cauchy problem for the scalar integro-differential equation 

/ poo \ 

rn ! ut(t,x)-[exp f(u x {t,y))dy) =0, u(0,x)=u{x). (1.1) 

O ■ \ Jx J x 

The model describes the slow erosion of granular flow, where u(t, x) is the height of the 
standing profile of granular matter. We assume that the slope of the profile has a fixed sign, 
i.e., u x > 0, and granular matter is poured at a constant rate from an uphill location outside 
the interval of interest, and slides down the hill as a very thin layer. The interaction between 
the two layers is controlled by the erosion function /, which denote the rate of the mass 
being eroded (or deposited if negative) per unit length travelled in x direction per unit mass 
passing through. We assume that the erosion rate / depends only on the slope u x . At the 
critical slope 1 (normalized), we have /(l) = 0. If u x > 1, we have erosion and / > 0. 
Otherwise, iiu x < 1, we have deposition and / < 0. The independent time- variable t denotes 
the amount of mass that has passed through, in a very long time. We will still refer to t as 
"time" throughout the paper, and call u{x) the "initial data". 

The model was first derived in pQ as the slow erosion limit of a 2 x 2 model for granular flow 
proposed by Hadeler & Kuttler in [TJ, with a specific erosion function f(u x ) = (u x — l)/u x . 
Later on, more general classes of erosion functions were studied, making distinction between 
whether the slope u x blows up or remains uniformly bounded. Let w=u x denote the slope. 
Under the following assumptions on the erosion function / € C 2 

/(1) = 0, /'>0, f"<0 (1.2) 
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the slope w remains uniformly bounded for all t > 0, see |2]. In this case, one could study 
the following conservation law for u>, 



Here the flux contains an integral term in x. Due to the nonlinearity of the function f(w), 
jumps in w could develop in finite time even for smooth initial data, which leads to kinks 
in the profile u. Thanks to the uniform bound on w, global existence and uniqueness of BV 
solutions for (|1.4|) are established in [21 [3]. 

However, if we allow more erosion for large slope, the solutions behave very differently. If 
the erosion function approaches a linear function for large w, i.e, if (j 1 . 3 1) is replaced by 



then the slope w could blow up, leading to vertical jumps in the profile, and w = u x would 
contain point masses. In this case one must study the equation (jl.ip . It is observed in [9 
that 3 types of singularities may occur in the solutions of (jl.ip , namely 

• a kink, where u x is discontinuous; 

• a jump, where u is discontinuous; 

• a hyper-kink, where u is continuous, but the right limit of u x is infinite, or both left 
and right limits of u x are infinite. 

The global existence of BV solutions for (jl.ip is obtained in [9], through a modified version 
of front tracking algorithm which generates piecewise affine approximations that also allow 
discontinuities. 

We remark that, model (jl.ip differs from other integro-differential equations in the lit- 
erature where the gradient u x may blow up. For example, the variational wave equation 
[H [8] and the Camassa-Holm equation [5] are both well-studied. In both cases, thanks to an 
a-priori bound on the L 2 norm of u x , the solution u remains Holder continuous at all times. 
In contrast, the solutions to our equation (jl.ip could develop jumps, and the distributional 
derivative u x could contain point masses. 

Since u(t,x) is an increasing function in x, the inverse X(t,u) is well-defined. We define 
the corresponding gradient 



Formal computation shows that X(t,u) and z(t,u) are conserved quantities, and they satisfy 
the equations 





w— >+oc 





z(t,u) = X u (t,u). 




(1.6) 



(1.7) 
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Here 

g(z) = zf(l/z) (1.8) 

is the erosion function in the coordinates (t, u), representing the rate of erosion per unit mass 
passed per unit distance in u. From the properties of / in (|1.2p and (jl.5|) . the function 
g(z) G C 2 satisfies 

5 (1) =0, g(0) >0, g"<0. (1.9) 

Note that, for a given t, when z(t, it) = on an interval in u, the physical slope w(t, x) 
blows up to infinity, and the profile u(t, x) has a vertical jump. However, the solutions for 
(|1.7| could become negative, which have no physical meaning. Therefore equation (|1.7p must 
be equipped with the pointwise constraint z > 0. We now modify equation (|1.7p into 



zt 



(g(z)-expj g(z(t,v))dv^j = M . (1.10) 



The measure in (jl.lOp yields the projection into the cone of non-negative functions. There- 
fore, we "transformed" the point mass in u x (t, x) into constraint in z(t,u). It turns out that 
the projection reduces the L 1 distance between solutions z(t,u). Thanks to this property, in 
[6] we proved continuous dependence on initial data and on the erosion function, for the en- 
tropy weak solutions generated as the limit of a front tracking approximation. This establishes 
a Lipschitz semigroup for the solutions of (|1.10p . 

In this paper we are interested in the traveling wave solutions of (jl.ip . We seek a traveling 
wave that connects profiles u(t,x) with slope w = 1 at both — oo and +oo, with slope w > 1 
in between, traveling with speed a, i.e., 

w(t,x) = W(£), € = x-<rt, lim W(£) = 1, > 1. (1.11) 

§— s-±oo 

In the variable u(t,x), these become profiles that travel upwards along lines of slope 1 with 
constant speed. See Figure [T] for an illustration. 

We now make an important observation. In Figured] we see that, viewed in the direction 
of the dotted lines with slope 1, the profile remains stationary in t. This motivates another 
variable change which would yield stationary traveling waves. 

To this end, we consider solutions z(t, u) to (jl.lOp such that, for all finite t > 0, 

z(t,-)£BV, z(t, O-leL 1 , lim z(t,u) = l, 0<z(t,u)<l. (1.12) 

u— >±oo 

This indicates that the profile u(t, x) approaches a linear asymptote with slope 1 at both 
x = dboo. Without loss of generality, we assume that u(t,x) approaches u = x at x — > +oo. 
We define the "drop" function 

r+oo 

q(t,u) = / (z(t, v) - 1) du. (1.13) 

J u 

Note that q(t, u) indicates the vertical drop of the profile u(t, x) at (t, u) comparing to the line 
u = x. Under our assumptions (j!.12[) . q(t,u) is an increasing function in u, and approaches 
as u — > +oo. Therefore, q(t,u) < 0. 

We also define the "total drop" of the profile as 

r+oo 

D = - / (z(t,v) -l)dv= \\z(t,-) - l|| L i. (1.14) 
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Figure 1: Example of a traveling wave profile for u(t,x). 



This denotes the vertical drop between the lines of asymptote at x — > ±00. See Figure [TJ We 
see that, as x — > —00, the profile u(t, x) approaches the asymptote u = x — D. 

We remark that z = 1 is a trivial solution. Under the assumptions (jl . 12[) . the L 1 norm of 
z — 1 remains constant in t. Therefore the total drop D also remains constant in t. 

If we assume further that z(t, u) = 1 outside an interval JcR, while z(t, u) < 1 on /, then 
the function u 1— > q(t,u) is invertible on I. Let u(t,q) be its inverse defined for q € (— -D,0). 
We now consider (i, q) as the independent variables, and define the composite functions 



In this new coordinate, the quantities X(t,q) and X q (t,q) are the conserved. We define the 
corresponding erosion function, 



X(t,q)=X(t,u(t,q)) 



C(t,q) =z(t,u(t,q)) . 



(1.15) 



h{z) = 



9(z) _ fO-M 



(0 < z < 1) 



Mi) 



</(!) = /'(I)- 



(1.16) 



l-z l/z-l 



For smooth solutions, we formally have 




(1.17) 



(1.18) 
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Treating C(t>l) as the unknown, and using the identity X q = — (), equation (|1.18p 
can be rewritten as 

(t-(l-0 2 (h(()-expj\(((t,s))ds^ = 0. (1.19) 

By this construction, smooth traveling wave solutions are stationary solutions of (|1.19p . 
This is the main technique in our analysis. We will construct traveling waves as stationary 
solutions of ((t,q) = Z(q) for (TTT91) . 

Depending on the size of the total drop D, different types of profiles can be constructed, 
and different types of singularities will form in the solutions. In particular, all the stationary 
profiles have a downward jump at q = —D (kink for the profile u), with possibly an interval 
where Z = (a shock in the profile u), and then possibly a smooth stationary rarefaction 
fan. In all cases, Z(q) is non-decreasing on — D < q < 0. See Section 3 for details. We will 
also show that such profiles are unique with respect to the total drop D. For any given D, 
there exist a unique traveling wave profile. 

The corresponding result for the physical variables u(t, x) and w(t, x) follows from the 
well-posedness of the variable changes (|1.15p . 

For the Cauchy problem of (jl.lOp with initial data satisfying (jl. 12[) . the existence of a 
Lipschitz semigroup of BV solutions is established in [6] . In turn this result provides us also 
the existence of semi-group solutions for the new variable £(£, q). 

We will study the local stability of the stationary traveling wave profiles. We show that, if 
the initial data is "non-decreasing", the solution approaches a traveling wave profile Z(q) as 
t +oo. By "non-decreasing", with a slight abuse of notation, we mean initial data which 
satisfies the following assumptions 

z(0,u) = < 1 ^ < Ua ' z (u 2 ) - % (ui) > for u 2 > u\ > u a . (1.20) 

[z {u) (u > u a ), 

Note that this property is shared by the traveling wave profile. As we will see later in Section 
4, property (|1.20p will be preserved in the solution for all t > 0. We will prove that, solutions 
of the Cauchy problem for (jl.lOp with initial data satisfying (I1.12P and (ll.20p . converges to 
the traveling wave profile asymptotically. Details will be explained in Section 4. 

The rest of the paper is organized as follows. In section 2 we make some basic analysis, 
where we derive waves speeds and prove some technical Lemmas. In section 3 we show 
the existence of traveling wave solutions by construction. Furthermore, such profiles are 
unique with respect to the total drop. We then return to the original variable and state 
the corresponding results in the physical variables u(t,x),w(t,x). In section 4 we establish 
local stability of these traveling waves, showing that solutions with non-decreasing initial 
data approach traveling waves asymptotically as t — > +oo. A numerical simulation is given 
in Section 5 to demonstrate the convergence. Finally, we give several concluding remarks in 
Section 6. 

2 Basic analysis 

2.1 Smooth stationary solutions for ((t,q) 

We start with the discussion on the properties of the erosion function h(z) defined in (I1.16p . 
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Lemma 2.1. For < z < 1, i/ie erosion function h(z) satisfies 

2h'(z) 

h(z)>0, h'(z)>0, h"(z)<^^-. (2.1) 

1 — z 

Proof. By the definition (|1.16|) . we have /i(0) = g(0) > 0. If < z < 1, since g(,z) > 0, we 
have /i(z) > 0. For z = 1, we have = — </(l) > 0. This proves that h(z) > and h(z) = 
if and only if when z = and <?(0) =0. 
For h' , we have 

T / / \ (1 ~ gV(g)+_g(g) 

( ' = (1 _ z )2 • ( 2 - 2 ) 

Since 

^{(l-z)g'(z) + g(z)} = (l-z)g"(z)<0, {(1 - *)</(*) + g{z)} z=l = 0, (2.3) 
we have /i'(z) > for < z < 1. For z = 1, we have 

ft'(l) = lim a-yW+PW = u «/(*) + *(*) = _ H y, {z) + h , {z)] = _ g , {1) _ m 

(1 — Z) z 1 ' 



so h'(l) = —^g"(l) > 0, proving the second inequality in ()2.1 
By £L2D and (|2"3j) . we now get 



^{(l-z)V(z)}<0. (2.4) 



Working out the derivative, we get 
d 



{(1 - z) 2 h'{z)} = -2(1 - z)ti(z) + (1 - z) 2 /i"(z) < . 



proving the third property in (I2.ip . □ 
For notation simplicity, in the rest of this paper we let F(£; q) denote the integral term 



F(C;q) = exp I h(C(t,s))ds, F t 
J a 



' q ((;q) = -h(()-F((;q). (2.5) 

Equation (|1 . 19[) can be rewritten as 

Ct - (1 - C) 2 (fc'(C)C - ^(0) • F(C; q) = 0. (2.6) 
For smooth solutions, along lines of characteristics 1 1->- q we have 

g(t) = -(1-cWO-nC;?), (2-7) 

C(t,9(*)) = -(l-() J fc 3 (().F(C; S ). (2.8) 

We observe that g(i) < and £(t,q(t)) < 0, therefore all characteristics travel to the left, and 
£ is decreasing along characteristics. 
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We now derive the ODE satisfied by smooth stationary solutions for (|2.6p . Let <fi(q) be a 
smooth stationary solution of (|2.6p . then it must be the solution of the Cauchy problem 



h 2 (6) 

The ODE (|2.9p is autonomous and can be solved explicitly by separation of variables. Indeed, 
we have 

h'(4>)4>' = h\4>), -+ ±h -> % = M- 

Integrating q over (q, 0), and /i over (h(<f>(q)), h (1)), we get 



^W" 1-^(1)9" 
This gives an explicit formula for 0(g), i.e, 



(2.11) 



where /i -1 denotes the inverse mapping of the function z i— > /i(z). By (|2.ip we know that 
h! > for < z < 1, therefore the inverse is well-defined. Furthermore, since 0' > 0, the 
function cj>(q) is strictly increasing. 

We observe that, if h(0) > 0, the solution (p(q) reaches at q = —D^k where is a 
finite value. By (|2.10p . we have 

We now define the function 

, f4>(q) for q £ [-D hk , 0] , 



(2.13) 

^0 for g < --D/jfe. 
We observe that if h(0) = 0, then D^k = +oo. If the total drop D is finite, we have 

> c , for q € ]-£>, 0] , where c D = <f>(-D) > . (2.14) 



Remark 2.2. If h(Q) = 0, this means we have g(0) = and /'(+oo) = 0. It is observed in [3| 
that the slope w = u x (t, x) remains uniformly bounded for all t. Correspondingly z(t, iz) ^ c Q 
for all (t, u) and for some positive constant c > 0. This implies 

C(*j <?) > c o > 0, Vt > and V<? € [-£>, 0] . 

We have the same observation in the traveling wave solution (|2.14p . 

We immediately have the next Lemma which provides a lower bound on the derivative of 
the smooth stationary profile. 
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Lemma 2.3. If h(0) > ; the smooth stationary profile <p satisfies 

h(Q) 2 

</>'(<?) >ci, Cl = V ; - - > 0, -D hk <q<0. (2.15) 

max 0<2 <i h'(z) 

If h(0) = 0, and let D be the total drop, we have 

<P'(q)>c 2 , c 2 = h{<t>{ ~ D)) \, ( - s > 0, -D<q<0. (2.16) 

Example 2.4. Let us choose the following erosion function 

f{w)=w--, g(z) = l-z 2 , h(() = l + C, 
w 

By ()2.1ip . the smooth profile satisfies 
We see that (£(-0.5) = 0, so D hk = 0.5. 

2.2 Wave speeds, admissible conditions and stationary waves 

To understand how the fronts move in the solution q), we derive here the wave speeds, for 
various types of singularities. Furthermore, we discus their admissibility conditions, following 
the Lax's entropy condition. We also single out the cases where the fronts are stationary. 
To simplify notation, we denote the integral term in the (t, u) coordinate by 

r+oo 

G(z;u)=ex.p g (z (t, v)) dv , G u (z; u) = -g{z) ■ G(z; u) . (2.17) 

J u 

Note that if C(^?) = - z (^ n ( ( z)) where q(t,u) is the drop function defined in (|1.13|) . we have 
F((;q) = G(z;u). 

Let u(t) be the location of a discontinuity, in a piecewise smooth solutions z(t,u). The 
wave speed ii(t) for three types of discontinuities in z(t,u) were derived in [6j. Let q(t) 
be the location of the corresponding discontinuity in the solution C(^<7)- Formally, by the 
conservation law (|1.7p we obtain 

q{t) = - [z(t, u+) - 1] u(t) - g (z(t, u+)) G (z; u) . (2.18) 

Thanks to (|2.18p . we can now list the corresponding wave speed q(t). 

Kink. A concave kink in u(t,x) corresponds to a downward jump in z(t,u). Since —g(z) is 
concave, only downward jumps are admissible. We can write 

z~ = z(t,u(t)-), z + = z(t,u(t)+), l>z~>z + >0. 

The speed of this wave is determined by the Rankine-Hugoniot condition for (|1.7p (see [6]), 

u{t) = -G{z-u)- 9{z ll-_ 9 z { f\ 
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By using the relation (|2.18p . we have 

9(f) = G (z; «) { (z + - 1) - g (z+ 

Using the functions h and F, we can write the speed as 

q(t) = -F (C; q) (1 - z+) (1 - z-) h ^_ = ^ +) . (2.19) 
Since z + < 1, the kink is stationary if and only if z~ = 1. 

Hyper-kink. A hyper-kink in w(t, x) corresponds to a downward jump in z(i, u) with z + = 
0. By taking the limit z + -> in (j2~T9l) . we get 

« t ) = -F(C;«)(l-,-)^^). (2.20) 

Note that a hyper-kink is admissible only if the jumps in z is downward. Furthermore, (|2.20p 
indicates that it is stationary if and only if z~ = 1. 

Shock. A jump in the profile u(t, x) corresponds to an interval where z(t, u) = 0. Let 
z(t, u) = on the interval u~ < u < u + , and we write 

z+(t) = z(t,u + (t)+) > 0, z~(t) = z(t,u~(t)-) > 0, A = u + -u~. 

Here A is the size of the shock. In the (t, q) coordinate, we let q~ < q < q + denote the 
corresponding interval of the shock, and A = u + — u~ = q + — q~ . This shock gives two fronts, 
namely the left front q~(t) and the right front q + (t). 
We now introduce the function 

h(o)s _ i 

ijj(8) = (»>0), V(0) = h(Q) > . (2.21) 

s 

The derivative of ifi satisfies 

^ = ^1 « ±1 > o, ( a > 0) (2.22) 

Let's first consider the right front q + (t). The speed u + (t) could be obtained by the 
Rankine-Hugoniot condition for (|1.1|) (see [BJ Subsection 2.1] for details) 

u + (t) = -^-±(g(z+)-iP(A)) . 



Using again (|2.18p . we get 



1-Z + 



Using the functions h and F, we get 



q + (t) = F (C; q + ) [V(A) - h (z + )] . (2.23) 
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A similar computation gives us the speed for the left front q , 

q -(t) = F (C; q~) ^(A)e- h (°) A - h (z")] . (2.24) 

We now discuss the admissible conditions. By Lax entropy condition, characteristics could 
only enter or stay parallel to shock curves. One can easily check that the left front at q~ is 
always admissible, see [6j [9] . For the right front q + , by using (|2.7|) . Lax condition yields the 
following inequality 

h(z + )-z + (l-z + )h'(z + )<il>(A). (2.25) 

Condition (|2.25p gives an upper bound for the value of z + for a given shock size A. To 
this end, it is convenient to define a function r/ (A) that gives the maximum value of z + such 
that the entropy condition f)2.25|) holds for a shock size of A, i.e., when fl2.25f) holds with 
equal sign. Therefore, we define the mapping A h > r] implicitly by 

h(rj) -rj)h!(rj) = ip(A). (2.26) 

We see that this implicit definition is well-defined. Indeed, we know tp' > 0. By the third 
property in (I2.ip we have 



j- (h(z) - z(l - z)h'(z)) = z [2h'(z) - (1 - z)ti'(z)] > 0, (0 < z < 1) . (2.27) 

We now check the condition when the fronts are stationary. The left front q~ is stationary 
if and only if z~ = 1, according to (|2,24p . For the right front q + , we define a function tp (A) 
such that for any given shock size A, we have q = if z + = 95(A). Thus, by the front speed 
(|2.23|) . the mapping A i->- (f is implicitly defined as 

h(ip) = %b(A). (2.28) 

Both h and ijj are strictly increasing functions, therefore the implicitly function (|2.28[) is 
well-defined. 

For a given shock size A, by (|2.23[) we have 

q + <0 if z + >ip(A), (2.29) 

q + > if z + < if (A) . (2.30) 

From (|2.23p we also observe that, if z + = 1, the right front of the shock is stationary, 
provided that it satisfies the admissible condition (|2.25p . Let D ss denote the smallest shock 
size such that z + = 1 is admissible. Then D ss satisfies the equation 

^{D ss ) = h{l) . (2.31) 

By (|2.22p , the function tp is strictly increasing, therefore the value D ss is uniquely determined 
by (I2T3TD . 

We can now conclude that any shock with z~ = z + = 1 and the shock size larger than 
D ss is stationary. 
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Example 2.5. If the erosion function has the simple form g(z) = 1 — z 2 , as in Example 12.41 
then h(z) = 1 + z and /i(0) = g(0) = 1. In this case, the implicit functions can be expressed 
explicitly. We have 

/e A - 1 e A - 1 e Dss - 1 

»?(A) = - 1 ' ^ (A) = — S~ _1 ' — 2^ — = 2 ' A. « 1-256. 

Let .D denote the total drop. For various values of -D, we plot the functions (j)(q), rj(D + q), 
(p(D + q) in Figure [21 using red, blue and green colors respectively. These graphs give us 
insights in the construction of the stationary traveling waves. Consider the intersection point 
of the graphs (f) (in red) and rj (in green), and let q be its q coordinate. Then, a shock 
with z = on q G [—D,q], with the right front value rj(D + q) would be stationary. From 
Figure [2] we see that, if D < D^, the two curves do not intersect, so no stationary shock 
could exist. If Dhk < D < D ss , then there exists only one intersection point, indicating one 
possible stationary shock. Finally, if D > D ss , the two curves do not intersect. However, 
since D > D ss , then a shock with z = on q £ [— D, 0] would be admissible with z + = 1, 
therefore it is stationary. 





o r ' 1 o.o H — 1 — 1 — 1 — 1 — • — 1 — ' — 1 — 1 — i"— ■ — 1 — 1 — 1 — 1 

-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 -0.2 -0.1 -0.0 -1.50 -1.00 -0.50 -0.00 



Figure 2: The graphs of <p{q) (in red), rj(D + q) (in blue) and ip{D + q) (in green) for 4 different 
values of D. (I): D < D hk . (II): D = D hk . (Ill): D hk <D< D ss . (IV): D > D ss . 

2.3 Technical Lemmas 

Inspired by the graphs in Figure [21 we now prove some technical lemmas. 
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Lemma 2.6. If h(0) > 0, the fallowings hold. 

a) The functions ip (A) and r\ (A) are strictly increasing for A > 0; 

b) ip(0)=v (0) = and <p (D ss ) = r) (D ss ) = 1; 

c) ip (A) < T) (A) for any A € (0, D ss ) . 

Proof. These properties follow from the definition of the functions (p and n, and the properties 
of the functions h and ip. 

a) Since h' > and ip' > 0, by the definition of ip in ()2.28j) we have cp' > 0. For r/ defined in 
(|2T26l) . by the property (|2371) and tp' > 0, we conclude that rf > 0. 

b) When A = 0, we have ^(0) = h(0). Since h(0) - (1 - 0) /i'(0) = h(0), from (1236]) we 
must have </j(0) = i](0) =0. 

When A = D ss , by d2~3TD and (^28]l . we have h(ip) = tp(D ss ) = h(l), therefore p(D ss ) = 1. 
For r/, by (|2.26|) we have h(n) = ip(D ss ) = h(l), therefore r/(D ss ) = 1. 

c) By the definitions (|2,26p and (|2.28p . and the property z(l — z)h'(z) > for < z < 1, we 
immediately conclude c). 

□ 

Lemma 2.7. If p(A) = <f)(q) for some values of A > and q S [—0^,0] then 

c/ ) '(q)-p'(A)> K >0, 

for some constant k > independent of A or q. This implies that any horizontal shift of the 
graph of p can intersect the graph of cp at most at one point. 

Proof. If /i(0) = 0, then (|2.28p implies h{p) = h(0), which gives p (A) = 0, so the Lemma is 
trivial, and one can let 

h 2 (z) 

k = mm . 

0<2<i h'(z) 

We now consider the case h(0) > 0. Let A and q be the values such that z = p(A) = (j)(q). 
By ([23]) and (I23HD . we have 

h'{pW{q) = h 2 (p) = V 2 (A), h'(p)p'(A) = f (A) . 



We have 



<P'(q)-p'(A) = -L[V, 2 (A)-V'(A)] > 3_^. mi? ^ (s) _^ (s)] 

h'(z) 1 J maxo< 2 <i h'(z) s>o 1 J 

1 . e h (°> (e h (°> - h{Q)s - 1) 



• mm 



maxo< 2 <i h'(z) s>o s 2 

1 e h <® 8 (e h( -°> - h(0)s - 1) 1 1 9 



lim i — — i = —— ■ -h 2 (0) > 0. 



maxo< 2 <i h'(z) s^-o+ s 2 maxo< 2 <i h'(z) 2 

We can now let 

h 2 (0) 

K= o ur-s >0, 

2 • max 0<z <i ft'(z) 

completing the proof. □ 
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Lemma 2.8. If h(0) > 0, we have Dhk < A 



Proof. By b) in Lemma [Z§] and P^jl . (gUP , we have 



D ss = 



Jo 




By Lemma 12771 [tp x ]' (£) > [0 x ]' (£), therefore Afc < As ; completing the proof. 



□ 



We immediately have the next Corollary regarding the intersection point of the graphs of 
(j)(q) and (p(q + D). This will be useful in the construction of the stationary traveling wave 
profiles and in the proof of their uniqueness w.r.t. total drop. 

Corollary 2.9. Let D be the total drop. We have the following results concerning the inter- 
section points of the graph of 4>(q) and (p(q + D) on the interval —D < q < 0. 

(1) . If D < Dhk, the two graphs never intersect, and (p(q + D) < <p(q) for —D < q < 0; 

(2) . If D = Dhk, the two graphs intersection once, at q = —D, and (p(q + D) < cj)(q) for 

—D < q < 0; 

(3) . If Dhk < D < D ss , the two graphs intersect once, at q + where —Dhk < Q + < 0, 

ip(q + D) > 4>{q) for q < q + and ip(q + D) < (j)(q) for q > q + ; 

(4) . If D = D ss , the two graphs intersect once atq = 0, and (p{q+D) > 4>{q) for —D < q < 0; 

(5) . If D > D ss , the two graphs never intersect, and (p(q + D) > cfi(q) for —D < q < 0; 
This Corollary is illustrated in Fig. [2j 

3 Existence and uniqueness of traveling waves 
3.1 Stationary traveling waves for ((t,q) 

In this subsection we prove the following Theorem on existence of traveling waves and their 
uniqueness with respect to total drop. 

Theorem 3.1. For every value of the total drop D, there exists exactly only one stationary 
traveling wave profile ((t,q) = Z(q), defined on q € [— D,0]. 

Proof. We prove the existence by construction. Let q) = Z{q) be defined on the interval 
q £ [-A0], with 



for various values of the total drop D. We recall the values Dhk and D ss , defined in (|2.12p 
and (|2.3ip . respectively. The stationary profile is constructed in different ways for different 
values of total drop D. 

Type 1. This is the case when D < Dhk- We let 



Z(-D) = 1, Z 



(0) = 1 



< Z(q) < 1 



q£]-D,0[, 




q = -D, 
-D < q < 0. 
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Here we have a kink at q = —D where Z has a downward jump, then it is connected to a 
smooth stationary profile on the right. By (|2.19p . the kink at q = —D is also stationary. 

We remark that if /i(0) = 0, the solution £(t, q) remains uniformly bounded away from 
for all t > 0. In this case, only Type 1 profiles are possible. For the next 3 types, we assume 
h(0) > 0. 

Type 2. In this case we have D = D^k- We let 



Z(q) 




q = -D, 
-D < q < 0. 



Here we have a hyper-kink at q = —D and it is connected to a smooth stationary profile on 
the right. By (|2.19p . the hyper-kink at q = —D is stationary. 

Type 3. We consider the case with D^k < D < D ss . Let (q + ,z + ) be the intersection 
point of the graphs <p{q) and <p(—D + q). According to Corollary 12.91 the two graphs intersect 
only once, at some interior point where —D < q + < 0. At this intersection point we have 

z + = (j)(q + ) = tp (A) , where A = q + + D . 

We define the profile 



Z(q) 



1 q=~D, 
—D < q < q^ 

[HQ) q + <q<0. 



Here we have a shock on q £ (—D,q + ), and it is connected to a smooth stationary profile on 
the right. The left front of the shock, located at q = —D, is stationary by (|2,24p . The right 
front at q + is also stationary by construction. 

Type 4. This is the last case when D > D ss . We let 

q=-D, 
Z(q) = { -D < q < 0, 

g = 0. 

This is a simple shock which is admissible. By (I2.24|) and (|2.23p . both left and right fronts 
are stationary. Therefore, the whole profile is stationary. 

Finally we note that the right front of a shock would be stationary if and only if it is an 
intersection point of the graph (j) with some horizontal shift of the graph (p. By Corollary 12.9^ 
there exist at most one such intersection point. This shows the uniqueness of these profiles 
with respect to the total drop, completing the proof. □ 

Example 3.2. If we choose the erosion function as in Example 12.41 then D ss satisfies 

2 = "P^"}- 1 , Dss „ 1.26. 
D 

^ss 

We have D^k = 0.5 from Example 12.41 The plots of these 4 types of traveling waves are given 
in Figure [3l 
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Type 1 



Type 2 




Figure 3: Stationary traveling wave profiles C(i, (?) = Z(q). 



3.2 Moving traveling waves for u(t, x) 

For the original physical variables u(t, x) and the slope w(t, x) = u x (t, x), the traveling waves 
are not stationary. They share some interesting properties which we would like to comment. 
Recall (jl.lip . and let W(£) be such a traveling wave, and let be the corresponding 

integral term. From the equation (jl.4p we get 

-^'K) + (»)^))' = o. 

Integrating it once, by the boundary condition (jl.lip we get 

- aW(0 + f(W(0)HO = -a, a= $™®\ ■ ^(0 (3.1) 
This gives us the constant propagation speed 

This says that, for any traveling wave that consists of a smooth part, (i.e., Type 1, 2 and 3), 
they must all travel with speed a = /'(l)! 

We now derive the ODE satisfied by W(£). Taking logarithm function on both sides of 
(PLED, we get 

ln<7 + ln(W- 1) =]nf{W) + J f(W(y))dy. 
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We now differentiate both sides of this equation with respect to £, and we get an autonomous 
ODE with a boundary condition 



w'(0 



f 2 (W)(W- 1) 

f(w)-f(w)(w-iy 



lim W(£) = 1. 



(3.3) 



Integrating once, we obtain the corresponding profile for u(t,x). 

If the shock size A is bigger than D ss , we have Type 4, which is a single shocks connected 
to u x = 1 on both the left and the right sides. It travels with the shock speed (see [HI (2.11)]): 



a 



3 A-/'(+oo) _ 1 

A ' 



As a consequence of Theorem 13.11 we conclude that such traveling waves exist for u(t, x) 
and w(t, x), and they are unique w.r.t. the total drop. 

Example 3.3. Consider the erosion function f(w) = w — ^ same as in Example 12.41 then 
([3.3p becomes 

w'(€) = -{w + if{w -i). 

Integrating W{£) in space, one obtains the traveling waves for the profile height u(t, x). These 
are plotted in Figure |4] for Type 1, 2 and 3, at t = 1. All these waves travel with the same 
speed a = /'(l) = 2. 



Type 1 



Type 2 



Type 3 




Figure 4: 3 Types of traveling wave profiles in the physical coordinate u(l,x) 
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4 Stability of traveling waves 



In this Section we prove the local stability for traveling waves. 

Theorem 4.1. Let z(t,u) be the solution of the Cauchy problem for (jl.lOp with initial data 
satisfying (|1.12p and (jl.20p . and let be defined in (|1.15p . Let D be the total drop, 

and Z{q) be the corresponding stationary traveling wave profile. Then, for every e > 0, there 
exists a finite value T e such that 



Here the constant C is independent of e. 

The rest of this Section is devoted to the proof of Theorem 14.11 

4.1 Structure of the solution ((t,q) for non-decreasing data 

Thanks to the additional assumption (|1 .20|) . we now make the following assumptions on the 
solutions q). 

Assumption 4.1. Let z(t,u) be the solution of the Cauchy problem for (| 1 . 1 U j) with initial 
data satisfying (|1.12j) and (|1.20p . Let ((t,q) be defined in (|1.15|) . It satisfies the following 
properties. 

• There is a stationary downward jump at q = —D, with —D-) = 1, for all t > 0; 

• There is maximum one shock in the solution, with the left front at q = —D, and the 
right front at q + {t) where —D < q + (t) < 0. 

• For any t > 0, q) is locally Lipschitz continuous and strictly increasing on the 
interval q(t) < q < 0, where q(t) = q + (t) if there is a shock, and q(t) = —D if no shock 
exists; 

• The total variation of (,{t,q) is uniformly bounded by 2 for all t > 0. 

Remark 4.2. The properties in Assumption 14. 1 1 are expected in the solution £(t, q). A rigorous 
proof could be carried out through front tracking piecewise constant approximate solutions. 
However, such a proof could be lengthy, since one has to enter the details of setting up the 
front tracking algorithm, then establish the a-priori estimates and the convergence. Since this 
current paper is focused on traveling waves and their properties, we would not provide the 
detailed proof and state these as assumptions. Below we provide some formal arguments to 
support these assumptions. 

(1). The assumptions fjl . 12f) and (ll.20p imply the following properties for C(0,<?), 



For smooth solutions £(t,q), the characteristic equations (I2.7p - ()2.8p hold. Therefore £ de- 
creases along characteristics. The assumption (jl.20p implies that ((t,q) < 1 for all t > and 
-D < q < 0. 

(2). The integral term F(£;q) is strictly decreasing in q, i.e, 



\\C(t,-)-Z(-)\\ hl <Ce, for all t>T e . 



(4.1) 



C(o, -D) = i, C(o,gi)<C(o,9 2 ) (gi<<a), C(o,o) = i. 



F(C;qi)> F((;q 2 ) 



(-D <q 1 <q 2 <0). 
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(3) . We consider two points in the solutions along the characteristics, £(t,qi(t)) and 
C(t,<fe(t)),with -D < gi(0) < q 2 (0) < and C (0,ft (0)) <_C(0,g 2 (0)). As long as ((t, qi (t)) > 
and C(t,Q2(t)) > 0, the equations (|2.7p - (|2.8|) hold. Let t be a first time such that 

0, qi(t)) = 0, q 2 (t)) or qi (i) = q 2 (t). 

Since the maps z i— > (1 — z) 2 h'(z), q i— > F (C; g) are decreasing and positive, we have 

J t [q2(t)-qi(t)]=q 2 (t)-qi(t)>0. 

This implies that (72 (i) = qi(t) cannot happen. If C(*>9i(*)) = C(*>?2(*)) holds, we can 
compute 

l[((i,q 2 (t))-at, qi (t))} =g(((t, qi (t)) [F(t; qi (t))-F((;q 2 (t))] > 

which gives a contradiction. In other words, for any t > 0, if ("(^ (/) > 0, then q 1— )• g) is 
non-decreasing and the rarefaction fan is strictly spreading. 

(4) . The total variation of •) is uniformly bounded by 2 for t > 0. 

(5) . If ((t,q) is strictly increasing, then there must be a downward jump in q) at 
q = -D, i.e., 

C(t,-D) = l, ((t,-D+)<l, (t>0). 

This jump is stationary. 

(6) . If /i(0) > 0, shocks might form in the solution. Since C(t,q) is strictly increasing, 
there could be maximum one shock, with the left front at q = —D and the right front at some 
q + where —D < q + (t) < for all t > 0. 

(7) . The solution q) remains smooth and strictly increasing on the part where ((t, q) > 
0. If there are no shocks, this is valid on the whole interval —D < q < 0. If there is a shock, 
then the interval is q + (t) < q < where q + (t) is the position of the right front of the shock. 

4.2 Properties of rarefaction fronts; formal arguments. 

Thanks to Assumption l4.lt we now only need to trace the location of the right front q + (t) of 
the possible shock, and study the evolution of the rarefaction fronts by characteristic equations 
(0)-(l2SI). Next Lemma follows from (l2T71)-(l2l 



Lemma 4.3. Let ((t,q(t)) > be a point on the rarefaction fan. Then, as t grows, the 
trajectory t 1— > (q, £) matches some horizontal shift of the graph of <j). Furthermore, the point 
(q(t),C(t,q(t))) travels to the left and downwards, until it merges into a singularity, at some 
t < T Co where T Co = {1 _ (o)mia %^ h , (z} , Co = C (0, q(0)) < 1. 

Proof. It suffices to observe that 

„ s h 2 (d>) C(t,q(t)) h 2 (() . . .. , „ 

therefore setting q Q = q(Q) < and Co = C (0) < 1) we have 

q(t) = - (1 - C (t, q(t))) h! (C (t, q{t))) F (C; g(t)) < - (1 - Co) mm < 0, 

0<z<l 

and the curve q{t) has to merge into a singularity before the time T^ o . □ 
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The formal arguments for asymptotic analysis is now rather simple, thanks to Lemma [4, 31 
We have the following observations. 

• Lemma 14.31 implies that, as t i— > +oo, the remaining rarefaction fan in the solution is 
generated near the point (q,z) = (0, 1). Again, since all rarefaction fronts travel along 
some horizontal shifts of <p(q), we see that the smooth part of the solution C(t>0) m ust 
approach the graph of 4>{q) as t — > +oo. 

• If a shock forms in the solution, by Corollary 12.91 and (|2.29|) - (|2.30j) . it must settle at the 
corresponding location of the stationary shock front in Z(q). 

Combining these two observations, one may conclude that £(t, q) converges to Z(q) in L 1 as 
t ->■ +oo. 

However, there is a complication. It is observed in [BJ [9] that, because of the admissible 
condition (|2,25p . characteristic curves could come out of the right front of the shock in a 
tangent direction at some t > 0. Therefore, the formal argument alone is not sufficient. 
Instead, in our proof we will construct suitable upper and lower envelopes for the solution 
C(t,q)- 



4.3 Upper envelopes. 

The upper envelopes might take 2 stages, depending on the type of stationary profiles. In 
stage 1 we control the smooth part of the solution. We show that the rarefaction fan gets 
very close to the stationary traveling wave after sufficiently long time. 

Lemma 4.4. Let D be the total drop, and let e > be given. We define the function 

'l q = ~D, 

+ (<l) = I 4>(q + e) —D < q < -e, 

1 -e<q<0. 

Then there exists a time Tf, such that for t > Tf , we have 

C(t,q) <4> + {q), -D<q<0. 

If the initial data satisfies £(0, q) < 0(g) for —e<q<0, we can simply take 4> + (q) = 4>(q)- 

Proof. Let e > be given and define q Q = — e < 0, Q Q = C,{Q,q ) < 1- We consider the 
rarefaction fronts generated on the interval (q a , 0). By Lemma [4 . 3 1 they travel along horizontal 
shifts of the graph of + , therefore they stay below this graph. Let t t— )■ q denote the 
characteristic curve with q(0) = q - By Lemma |4~3| the point (q(t),C(t,q(t))) will merge into 
a singularity before the time Tf = T^ o . After that, the smooth part of the solution is the 
rarefaction fan generated on the interval (q o ,0) at t = 0. Therefore, we have C(t,q) < 
for t>T{. 

Finally, if £(0,9) < 4>(q) for — e < q < 0, we can simply take (p + (q) = <p{q) and carry out 
the whole argument in a complete similar way, completing the proof. □ 



If the stationary profile is of Type 1 and Type 2, the upper envelopes are complete. We 
now consider Type 3 and Type 4, and enter Stage 2 to control the location of the right front 
of the shock. 
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Lemma 4.5. Assume D > Dhk, and 1st q + < be the location of the right front of the shock 
in the stationary traveling wave profile Z{q). Let e > be given. There exists a time Tf , such 
that 

C(t,q)<Z+(q), -D<q<0, t > Tf + Tf , 
where the function Z + is defined as 



Z+(q) 



1, q=~D, 

0, -D<q<q, q = q +-C 1 e, (4.2) 



> + (q), q<q<0- 
Here the constant C\ does not dependent on e. 

Proof. Let D > Dhk be the total drop, and let q be the intersection point of the graphs of 
4> + (q) and ip(q + D), such that <j) + {q) = f(q + D). By Lemma 12.71 the functions (j) and (p are 
strictly increasing and transversal. We have 

+ - - max 9 <//(<?) 1 h 2 (z) 
q + ~q<Ce, C = H - = - max — — , 

K K 0<z<l h'(Z) 

where k is defined in Lemma |2,71 We can choose the constant in (|4.2[) as C\ = 2C. 
We now construct the upper envelopes, for t >T[, 



Z + {t,q) 



1, q = ~D, 

0, -D< q < q(t), 

l>+(g), q(t)<q<0, 



Here the front q{t) travels with the speed as if it were the right front of a shock. By (|2.23p . 
we have 

q'(t) = F(C; q) ■ 1 ~ + ( "t } (g) [^(q + D) - h(<f>+(q))] , q(T() = -D hk - e . (4.3) 

The ODE (14. 3D is solved for t > Tf, until at some time t = Tf + Tf when the front q(t) passes 
the one in Z + , i.e., when q(Tf + T| ) > q + — C\E. 

We now show that T| is finite for any given e. When q(t) < q the graph of cj) + lies strictly 
below the graph of ip(q + D). Since q + — C\£ < q, therefore by continuity we have 

v e = min {q + D) — cf> + (q) : q < q + — C\e} > 0. 

Hence, as long as q(t) < q + — C\£, we can compute 

q\t) > [l-<f> + (q)] [^(q + D)-^^^))} 

> [l-Hq + )] [h{ip{q + D))-h{ct> + {q))] 

> [l — </ ) (q ,+ )] • min/i' • v £ . 

We can now conclude that 

D hk + e 



Z + (t, q) < Z + (q), for t>Tf+ Tf , where Tf 



[1 — <p(q + )] ■ min/i' • v £ 
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It remains to show that the solution C(t,q) satisfies 

at,q)<Z+(t,q) for Tf < t < Tf + If . 
Indeed, by construction (|4.4|) holds for t = Tf because 

Z+(Tf,q) =<P + (q)>((T(, q ), 



(4.4) 



where <fr is the the function defined in Lemma [4.4i Now we consider a later time i > Tf. By 
Lemma |4.3| the smooth part of the solution <^(t,q) remains below the graph of (f) + {q). It is 
then enough to check the speed of right front of the possible shock. Assume that C has 
a shock, with its right front at q(t), and 

q(t) = q(t), C(t,q(t}+) < </> + (^~))- 

By ()2.23p we clearly have q'(t) > q'(t), so the graph of C remains below that of Z + (t,q), 
completing the proof. □ 

The two stages are illustrated in Figure [5] for Type 3. 




Figure 5: Upper envelope, for type 3, stage 1 (left) and stage 2 (right). 



4.4 Lower envelopes. 

For Type 4, the lower envelope is trivial, since one can simply take the stationary traveling 
wave profile Z{q). We now construct the lower envelopes for Type 1, 2, and 3, following 
a similar line of arguments as those for the upper envelopes. In stage 1, we show that the 
rarefaction fan approaches the stationary traveling wave profile as t grows. 

Lemma 4.6. Assume D < D ss and let e > be given. We define the function 



q = -D, 
—D < q < qi, 
qi<q<0, 



(4.5) 



Here the value q\ is determined as follows. If D is sufficiently small such that the graph of 
<j)(q — e) lies completely above the graph of w(q + D), we let q\ = —D, and we remove the 
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second line with —D < q < qi in the definition (|4.5p . Otherwise, we let q\ be the right-most 
intersection point of those two graphs. Then, there exists a time T-[ , such that, 

C(t,q) >4>-{q), t>T{. (4.6) 

Proof. The proof follows a similar structure as the one for Lemma 14.41 with some modifica- 
tions. Let e > be given. Since 4>~(0) = <f>(— e) < 1, by continuity there exists q Q < such 
that 

C(0,<Z) ><T(g), for qo<q<0, ( a = C (0, q a ) < 1. 

We consider the rarefaction fronts generated on [q a , 0] and let t h > g be the characteristic 
curve with g(0) = q . Then, the point (q(t)X(t,q(t))) travels along some shift of the graph 
of <j> until it merges into a singularity. There are two possibilities. 

(1) . If D is small such that the graph of <f> lies completely above the graph of r](q + D), 
then the characteristic will reach q = —D and enter the downward jump at q = —D. 

(2) . If the two graphs intersect, then the right-most position of the right front of a possible 
shock is qi, and (|4.6|) holds. If the actual shock front is to the left of q%, or the characteristic 
does not enter any shock, we still have (|4.6p . 

By Lemma [4. 3 1 we have T\ = T^ , which is finite because C, Q < 1, completing the proof. □ 



If D is small such that q\ = —D in Lemma l4.6[ the lower envelopes are complete. For the 
rest of the subsection, we assume q\ > —D. We now enter stage 2, and we control the location 
of the right front of the shock. We will use again the symbols q, q, q, C etc, for different values, 
without causing any confusion. 



Lemma 4.7. Assume D < D ss and q\ > —D in (J43J). Let q + be the location of right front 
of the stationary shock for Type 3, and let q + = —D if it is Type 1 and 2. Let e > be given. 
There exists a time 7J , such that 

t(t,q)>Z-(q), -D<q<Q, for t>T{ + T%. 

Here the function Z~ is defined as follows. For Type 1, we let 

Z-{q) = 

For Type 2 and 3, we let 



1, q = ~D, 

-e), -D<q<0. 



Z-(q) 



1, q = ~D, 

0, -D<q<q, q = q + + C 2 e, (4.7) 



>(g-e), q<q<0. 
The constant C 2 does not depend on e. 

Proof. Again, the proof follows a similar line of arguments as for Lemma 14.51 with modifica- 
tions. Let q be the intersection point of the graphs of <p{q — e) and (p{q + D). We have 

- 1 h 2 (z) 
q — q + < Ce, C = — • max 



k o<z<i h'(z) 
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We can now choose the constant in (|4.7p as C 2 = 2C. 
The lower envelopes are defined as follows, for t >Tf, 



Z~(t,q) 



1, q=~D, 
0, —D < q < q(t), 

{(j>(q-e), q(t)<Q<0- 



Here the front q(t) travels with the speed as if it were the right front of a shock. By (|2.23p . 
q(t) satisfies the ODE, for t > T{ 

q'{t) = F((; q) ■ 1 ~ ^ ~ g) [^(q + D) - h(</>(q - e))} , q{Tf) = qi . (4.8) 



The ODE (|4.8p is solved for t > Tf, until at some time t = T{ + T 2 when the front q(t) passes 
the one in Z~ , i.e., when q{T{ + T{) < q- 

We now show that 7J is finite for any given e. As in the proof of Lemma 14. 5\ observe 
that when q(t) > q the graph of <f>~~ lies strictly above the graph of (p(q + D), therefore, by 
continuity and since q~ + C^e > q we have 

v £ = — max [if (q + D) - (f>~(q) : q > q + + C2e} > 0. 

Hence, as long as q + + C*2£ < < qi, we can compute 

< [l - [^(g + £>) - 
<[l-0(?i)] + 

< — [1 — <p(qi)] • vahih • v £ . 

We can now conclude that 

D hk + e 



Z~(t, q) > Z~(q), for t>T{ + T 2 £ , where T 2 £ 



[1 — 0(^i)] • min/i' • u e 



We note that for Type 1, the front q(t) would merge into q = —D at some t < Tf + T 2 ■ 
It remains to show that the solution £(t, q) satisfies 

C(t,q)>Z-(t,q) for Tf < t < T( + T 2 S . (4.9) 

Indeed, by construction (14.90 holds for t = T{, since Z~(T{,q) = <f>~(q), where <p~ is the the 
function defined in Lemma 14.61 Now we consider a later time i > T{ ■ By Lemma 14.31 the 
smooth part of the solution q) remains above the graph of (j){q — e). It is then enough to 
check the speed of right front of the possible shock. Assume that z(t, q) has a shock, with its 
right front at q(t), and 

q(t) = q(t), ((i, q(t)) > (j)(q(t) - e). 

By (|2.23p we clearly have q'it) < q'(t), so the graph of C(*><z) remains above that of Z~(t,q), 
completing the proof. □ 

The two stages are illustrated in Figure for Type 3. 
We now combine all the estimates and prove Theorem 14.11 
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Figure 6: Lower envelope, for Type 3, with Stage 1 (left) and Stage 2 (right). 



Proof, (of Theorem 14. ip Given e > 0. By Lemma 14.51 and Lemma 14.71 there exists a time 

T £ = max{Tf + T| , 7? + 7?} , 

such that for t > T £ , we have 

< C(t, q) < Z+(q), -D < q < . (4.10) 
Here Z~ and Z + are defined in (j4.7|) and Q4.2j> . respectively, and they satisfy 

< < Z+(q), -D<q<0, (4.11) 



^IIti <c l£ , 



Z-|| T1 <(Ci + C 2 )e, (4.12) 



where the constants Ci, C 2 does not depend of s. Thanks to (|4.10p - (|4.12p . we now conclude 

||C(t, •) - Z(.)|| L i < \\Z+ - Z"|| L1 < (Cx + C 2 )e, t > T £ , 
completing the proof. □ 



5 A Numerical Example 



In this section we present a numerical simulation of (jl.lOp . We generate piecewise constant 
approximate solutions using an extended version of the wave front tracking algorithm de- 
scribed in [6j [9] . We use the erosion function 

1 



g(z) = (I -z) [-+:. 



and the initial data 



z (u) 



1 


lexp 



-\{u 



for u < , 
for < u < 0.6 , 
+ 0.11)) for 0.6 <u. 
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With this initial data, the traveling wave profile is of Type 3. Solutions are plotted for 
nine different values of t, for both the functions C(*>9) an d z(t,u), in Figure [7] and Figure 
respectively. In Figure [7] we also plotted the stationary Type 3 traveling wave Z(q) (in red) 
together with the solution ((t,q). 

As we observe in Figure El the traveling waves for the solutions z(t,u) are not stationary. 
We clearly observe the moving traveling wave in the last 2-3 plots in Figure It is interesting 
to observe that, for waves of Type 1, 2 and 3, they all travel with the same constant speed a 
which can be deduced from (13.21). i.e., 



1 

a 



/'(I) 



3 
2 



1.5. 



For a Type 4 wave, it travels with the shock speed. 

This simulation also demonstrates the complexity of the transient dynamics of the wave 
formation and interaction. One obverses that the shock in the initial data disappears as the 
rarefaction wave on the right merges into the shock, only to reappears later as a new shock 
forms. 




Figure 7: Plots for the function ((t, q) vs q for various values of t. 
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t = 0.0000 



t = 0.0842 






Figure 8: Plots for the function z(t, u) vs u for various values of t. 



6 Concluding remarks 

In this paper we prove the existence of traveling wave profiles for an integro differential 
equation modeling slow erosion of granular flow. Such profiles are unique with respect to 
the total drop D. Furthermore we show that these profiles provide local attractors for the 
solutions of the Cauchy problem. 

We now conclude the paper with several final remarks. 

Remark 6.1. The basin of attraction of the traveling wave profile is actually much larger, 
and the initial data does not need to be non-decreasing. The initial data Co(?) = C(0, q) only 
needs to satisfy the following. For some e > 0, we have 



Co(--D) = Co(0) = 1, Co(9) < 1 - Ce, 

Cote) - Cote) > o, 
Cote) - Cote) > o, 

TV{Co}<M, ||Co - l|| L i < M. 



-D + e<q<-e) 
-e<q 2 <qi< 0) 
-D < g 4 < q 3 < -D + e) 



(6.1) 
(6.2) 
(6.3) 
(6.4) 



From [6] we see that for general initial data with bounded variation, the total variation 
of C(ij ■) can grow exponentially in t. However for this simpler case (|6.ip - (|6.4p . one should be 
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able to improve the BV estimate and actually obtain a bound that is uniform in t. We now 
provide a formal argument. Consider a characteristic curve t i— > q(t) initiated on the interval 
[-D + e, — e]. As long as ((t,q(t)) > 0, we have q < —cqc and C(t,q(t)) < —cq€, i.e., the 
characteristic curve travels strictly to the left, and the £ value is strictly decreasing along the 
characteristics. In finite time, this curve will either enter a shock such that £ = 0, or reach 
q = —D. This implies that all the singularities would finish all possible interactions in finite 
time, and after that the solution C(^>(?) wm De non-decreasing, as in the assumption (|l,20p . 
Therefore, after finite time the total variation of £(t, •) will be bounded by 2. Then one can 
apply the result in this paper and obtain the asymptotic behavior. 

Remark 6.2. At this point, we also conjecture that our result could be extended to general 
BV initial data z (u), provided that the total drop D is positive. Due to the nonlinearity of 
the erosion function, all singularities will interact and merge into a single singularity in finite 
time, even though the transient dynamics could be very complicated. The solution will satisfy 
the assumption (|1.20p in finite time. It should be possible to carry out a rigorous analysis 
through piecewise constant approximate solutions generated by the front tracking algorithm. 
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